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Abstract 

We consider the thermodynamics of a black hole coupled to thermal 
radiation in a spatially finite (spherical) region. Thermodynamic state 
functions are derived in the canonical ensemble, defined by elements of 
radius Tq and boundary temperature T(ro). Using recent solutions of 
the semi-classical back reaction problem, we compute the 0(?i) correc- 
tions to the mass of the black hole, thermal energy, the entropy and free 
energy due to the presence of hot conformal scalars, massless spinors 
and U(l) gauge quantum fields in the vicinity of the hole. The free 
energy is particularly important for assessing under what conditions 
the nucleation of black holes from hot flat space is likely to occur. 
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1 Introduction 



A Schwarzschild black hole in empty space radiates quanta possessing a tem- 
perature characterized by the mass M of the hole. At large distances from 
the hole (r ^ M), the temperature of the radiation approaches Too = g;^- 
Studies seeking to relate the thermal properties of black holes to quantum 
gravity typically omit the radiation from consideration. Nevertheless, it is 
clear that black holes are neither thermally nor mechanically isolated from 
their surroundings and in deducing the thermodynamic properties of black 
holes, the presence of the radiation field should be taken into account. One 
way in which to do so is to solve the semi-classical Einstein equation 



taking as source term the gravitationally induced renormalized stress ten- 
sor < Tf^i, > characterizing the thermal radiation. The calculation of the 
modified spacetime metric is the back reaction problem, and it is expected 
that the new metric gives a better approximation to the spacetime geome- 
try associated with thermal equilibrium than one satisfying the source-free 
(< T^u >= 0) Einstein equation. In equilibrium, the ensemble must be time 
independent, thus the perturbed metric must be static, as well as spheri- 
cally symmetric, as the stress tensors employed here are renormalized on a 
Schwarzschild background. From knowledge of the modified metric, one can 
calculate the 0{h) correction to the black hole temperature [|^, the entropy 

and the effective potential [^]. In the present paper, we use solutions of 
the back reaction problem obtained in Q, Q and in order to compute 
the fractional corrections to the mass, the thermodynamic energy, the free 
energy and the entropy for a black hole in the canonical ensemble arising 
from its interaction with scalar, spinor and U(l) gauge quantum fields. 

An important feature in this analysis is the use of finite spatial bound- 
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aries. Indeed, it was recognized in [T| that the back reaction problem (1) 
has no definite solution unless the system composed of black hole plus radi- 
ation is enclosed in a spatially bounded region or "box" . This is due to the 
fact that the stress tensors renormalized on a black hole background space- 
time are asymptotically constant, so the corrections to the metric do not 
remain uniformly small for sufficiently large radius. Physically, this means 
the radiation in a box that is too large would collapse onto the hole thereby 
producing a larger one. Yet, there is considerably more to the introduction 
of finite spatial boundaries than just as a means for ensuring perturbative 
validity of the solution. As pointed out in there are additional advan- 
tages to be gained by employing boundary conditions at a spatially finite 
location, as opposed to spatial infinity. In the first instance, as already 
noted, self-gravitating thermodynamic systems are not usually asymptoti- 
cally flat. In particular, in the canonical ensemble, the system boundary is 
to be maintained at a constant temperature, and this is achieved by coupling 
the boundary to an external heat reservoir. Such an arrangement obviously 
fails to satisfy asymptotic flatness (for the region exterior to the system pos- 
sesses a constant energy density), but with finite spatial boundaries there is 
no need to assume asymptotic flatness in the spatial directions. Secondly, 
the usual thermodynamic limit requiring infinite spatial extent does not exist 
for equilibrium self-gravitating systems at finite temperature. This follows 
since the system is unstable to gravitational collapse, or recollapse if a black 
hole is already present. This in practice presents no problem since physi- 
cally, one only requires that the system can in principle become sufficiently 
large so that the fluctuations become negligible. 

This instability is reflected in the formally negative heat capacity, which 
in turn, implies a divergent canonical partition function 1^]. Enclosing the 
black hole in a box has the effect of stabilizing the hole yielding a system with 
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a positive heat capacity and allows one to derive a meaningful, convergent 
partition function . The key to these results is the fact that the black hole 
energy and temperature at a finite radius are not inversely related, as they 
are at spatial infinity, because of the blueshifting of temperature in a static 
gravitational field. 

Given the advantages of invoking finite spatial boundaries, the choice 
must be made in specifying the type of boundary data which is to be fixed. 
In the present work, we consider the canonical emsemble, in which the tem- 
perature is held fixed at the cavity wall. The elements of this ensemble thus 
consist of spherical cavities of radius ro and temperature T{ro) enclosing a 
black hole at their centers and filled with thermal radiation. The boundary 
temperature is held fixed by coupling the cavity wall to a large external heat 
reservoir. At first sight, this might seem like an artificial arrangement, how- 
ever the heat reservoir is intended to represent the natural universe and the 
boundary conditions mimic the implantation of the hole and its equilibrating 
radiation into the universe. 

One of the main objectives in this paper is the calculation of a thermody- 
namic potential appropriate for a black hole coupled to a thermal radiation 
field. This motivated by the fact that the nucleation of a black hole from 
hot flat space (gravitons and other massless fields on a flat background ge- 
ometry) can be viewed as a phase transition In the canonical ensemble 
the temperature and volume (area of the cavity wall) are fixed, so the rele- 
vant potential is the Helmholtz free energy F = E — TS, E is the thermal 
energy and S is the entropy. This is related to the canonical partition func- 
tion by I3F = — log(Z), j3 = T^^. In a phase transition, the value of F 
should decrease and if the free energy is a (local) minimum, the system is in 
a state of (meta) stable equilibrium. To ascertain the likelyhood of nucle- 
ation, one compares the free energies evaluated for hot flat space and black 
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hole configurations. Such an analysis has been addressed in for a free 
black hole. 

We begin in Section II by deriving the mass corrections for a black hole in 
a finite cavity of radius Tq valid to 0(e) in the back reaction, where e ~ (^)^j 
and Ipi is the Planck length. The action / of the equilibrium metric is 
calculated in Section III and depends on the trace anomalies of the stress 
tensors used in (1). The thermodynamic energy, entropy and the free energy 
of the combined system comprised of black hole plus radiation are derived 
assuming that the partition function contains / as its leading term. The 
0(e) fractional corrections to the energy and entropy are calculated to reveal 
the spin dependence of the back reaction. In Section IV we discuss some 
aspects of black hole nucleation based on the results of these calculations. 
The explicit solutions of the back reaction (1) used here are collected in an 
Appendix. 



2 Schwarzschild Mass in a Spherical Cavity 

It is well-known that the equilibrium temperature distribution of a static, 
self-gravitating system is not constant, but is blueshifted, indicating that 
temperature is actually scale-dependent and is thus, not a purely intensive 
thermodynamic variable. This distribution for a Schwarzschild black hole 
including quantum mechanical back reaction is given by |||],[0] (li; = ^^) 



T{r) = (8^M)-i(l-u;)-^/2 



1 - e{p{w) - riK-^ - ^(1 - wy^niw)} 



(2) 

(see the Appendix for the definitions of p, n and p and K) where e = 
(^)2 < 1, Mpi = h^/"^ is the Planck mass (in units where G = c = kp = 1). 
An important finding is that the Schwarzschild mass is double- valued in the 
canonical emsemble O, p. Indeed, since the local temperature is fixed at 
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the boundary radius of the spherical cavity, we can invert (2) and solve for 
M = M{ro,T{ro)). By squaring (2) and remembering that (we henceforth 
drop the subscript "o" from the remainder of this discussion) 

le^MI ^ \e{p - nK-' " f (1 " ^)"V}| < 1, (3) 

we obtain the equation, valid to 0(e) 

- + = 2ecj2j^(tt;), (4) 

where a = ^^^j. > 0. Although this is a nontrivial transcendental equation 
in w {J- contains positive and negative powers of w and a logarithm), we 
can solve for approximate solutions by perturbing around the exact solution 
which has been computed for e = 0. In the absence of any back reaction, it is 
known Q that there are no positive solutions when rT < \/27/8vr ~ 0.207. 
When rT > ^/27 / Bvr there are two real nonnegative solutions of the form 

wi = ' =-[l-2cos(- + -)], (5) 

2M2(r,T) 1, rt,, 
W2 = ^' ' =-[l + 2cos(-)], (6) 

COs(q) = 1 - yCJ^ (7) 

< a < vr, (8) 

with M2 > Ml. For rT = ^/27/87r (a = vr), 3M2 = 3Mi = r indicating 
the cavity wall coincides with the circular photon orbit, while for rT 00 
(a — > 0) we have M2 § and Mi — > ( see Fig. 1 ). For the heavier 
mass branch (M2), the cavity radius can run from the black hole horizon 
out to the circular photon orbit. The mass increases as the wall is pulled 
in. For the lighter mass. Mi decreases as rT — > 00. By virtue of the 
parametrization of the solutions (5) and (6) in terms of the angle a, the 
entire semi-infinite two-dimensional parameter domain \/27 /8vr < rT < 00 
of the r — T plane can be compactified to a finite one-dimensional interval. 
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This makes it particularly easy to see, at a glance, how changes in the 
cavity radius and wall temperature affect the value of the black hole mass 
consistent with the given boundary conditions. Thus, for given r > and 
r > such that rT > \/27/87r, one calculates a through (7) and then the 
values of Mi or M2 can be read off from the respective branch curves in 
Fig. 1. Moreover, this plot also indicates the allowed range in the variable 
w = 2M/r, which appears in all our back-reaction formulae. For example, 
for the heavier mass solution {AI2) wc see immediately that 2/3 < W2 < 1, 
whereas < wi < 2/3 for the Mi branch. This is important, for it means 
when we come to evaluate some function of w, we should restrict the domain 
of that function to correspond to the range of w. This summarizes the 
thermal and radial dependence of the mass of the vacuum (that is , without 
the back reaction) black hole in the canonical ensemble. We will see shortly 
how the back-reaction modifies the allowed range of wi^2- 

We now couple the hole to the ambient thermal radiation in the cavity 
and ask for the fractional corrections to the mass of the hole arising from the 
back reaction. If (3) holds over some range in w, we can obtain approximate 
solutions to (4) by assuming an ansatz of the form 

Siv 

w = w + Sw, I — 1<1 (9) 
w 

for the perturbed solutions. Inserting (9) into (4) and retaining only the 
0(e) and 0{Sw) terms yields 

fa= Y^W. (10) 

w{3w - 2) ^ ^ 

The approximate solutions of (4) expressed in terms of the angle a are thus 

~ f ^ ( ^ , ^€1,2(1 -cos(a))J^(wi,2(«)) 

u;i,2(Q;)(3u;i,2(a) - 2) 

where we have used (7) and (9) and wi^2 refers to one of the solutions (5) or 

(6). In determining the size of the mass perturbations, it should be noted 



that e itself depends on the angle a: 




Mpi 




) 



2 



(12) 



M2 



|[l + 2cos(f)] 



(with a similar expression for ei). Thus, in varying the product rT or a, we 
can hold e < 1 to a fixed value by varying r accordingly. So, for example, 
we can maintain M2 constant as rT — > co by simultaneously decreasing r 
as indicated in Fig. 1. This corresponds to a curve of constant mass in the 
r — T plane. On the other hand, if r is fixed, then varying a corresponds 
to varying the temperature at the cavity wall. In any case, we can ensure 
that e < 1 by taking the cavity radius to be r > 31 pi for the heavier mass 
solution. The pure perturbation term Aw = Jf^^J^2) fro™ (H) is plotted in 
Fig. 2 and Fig. 3 for the M2 and Mi mass branches, respectively. There 
is a striking qualitative as well as quantitative contrast among the three 
types of back reaction we are considering here. First, it should be noted 
that in all cases there is a divergence in Aw as a — > vr corresponding to the 
wall radius approaching the circular photon orbit (coincidentally, the specific 
heat at constant wall area also diverges at r = 3M, both with and without 
taking the back reaction into account). This pole is manifest in (11) and 
simply means that the mass perturbation cannot be extrapolated reliably to 
this limit. Without having the exact solution(s) of (4), it is difficult to judge 
whether this pole is a significant feature of the back reaction or is merely 
an artifact of the approximation. We can only say that reliable results will 
require calculations beyond 0{h). Nevertheless, there is a wide range in a 
for which the perturbation term Aw is truly small and hence consistent with 
the expansion (9). This range depends, as indicated by the Figures 2 and 
3, on the spin of the quantum field. While the spinor tends to increase the 
mass of M2 as a ^ vr, the scalar and gauge particles tend to lower it. The 
perturbations valid for the lower mass branch behave in quite the opposite 
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manner. That is, while the spinor back reaction now tends to lower the 
mass Ml from its vacuum value, the scalar and vector boson tend to increase 
it. Apart from the pole at wi^2 = 2/3, all three curves for the Mi branch 
diverge as a — > 0. However, this divergence is well understood, and is due to 
the fact that the functions characterizing the metric perturbations (/? and /u) 
grow without bound as r — > oo |l|, ^, ^, ^. This in turn is a consequence of 
the fact that the stress tensors renormalized on a Schwarzschild background 
are asymptotically constant, as stated previously in the Introduction. So, 
in this small a regime, a natural cut-off is provided by the cavity radius. 
Referring back to the M2 branch plus perturbations, we see that the lower 
limit in the range in W2 is extended towards by the effects of the scalar 
and gauge boson back reactions. Clearly, as a — > vr, these perturbations 
become arbitrarily large and negative and would result in an M2 < 0, which 
is unphysical, as well as being outside the domain of perturbative validity 
(for \Aw\ > 1). The spinor tends to extend the upper limit in the range of W2 
to be greater than 1, but this would translate into a cavity radius inside the 
black hole event horizon. Thus, this perturbation must be cut off before a 
reaches vr. In the case of the lower mass branch, provided we remain within 
the range of perturbative validity, we note that the spinor perturbation tends 
to lower this mass while the scalar and vector tend to increase it. Except 
for the behavior close to a = vr, the conformal scalar field contributes a 
negligible amount to the mass corrections for both branches in comparison 
to the spinor and gauge boson. Moreover, these latter two fields give rise 
to competing effects, as the signs of their associated mass perturbations are 
opposite. This sign difference becomes all the more important when effects 
of multiple-field back reaction are treated. Multiple particle species arise, 
for example, in gauge theories of particle physics. The actual number of 
(fundamental) scalars, spinors and gauge bosons depends on the symmetry 
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group and the dimensionality of the group representations. 

Perhaps the most important point to be drawn from all this is that it 
is possible for the roles of M2 and Mi to be dynamically switched, that is, 
although at leading order we have M2 > Mi , the effects of the back reaction 
may result in a level crossing such that Mi > M2. Thermodynamically, 
it is the heavier of the two masses which gets nucleated so the back 
reaction must provide an important ingredient in assessing the likelyhood 
of the nucleation of black holes from hot flat space. 

3 Action, Thermal Energy and Entropy 

In deriving gravitational thermodynamics from a Euclidean path integral, 

Z = j d^l[g,ct>\e-'^3M, (13) 

one expects the dominant contribution to the canonical partition function to 
come from metrics g and fields (j) that are near a background metric g^^^ and 
background fields <t)^'^\ respectively. These background fields are solutions 
of classical field equations. First-order quantum effects can be incorporated 
in Z by expanding the action in a Taylor series about the background fields 

g = gi^) + ~g^ <^ = ^(0) + ^, (14) 

SO that 

I[g, 0] = I[g^^^] + l2[g\ + l2[4>] + higher order terms, (15) 

where I2 is quadratic in the field fiuctuations, and we have set the back- 
ground matter (or radiation) fields to zero, this corresponding to the case at 
hand (i.e., Hawking radiation is quantum mechanical, not classical). As is 
well known, the functional integration of I2 with respect to g and (j) leads to 
determinants of differential operators which can be exponentiated formally 
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to yield, when added to /[f?'-'^^], the one-loop effective action. However, from 
the point-of-view of the back reaction of quantum fields on the background 
geometry g^^\ the solutions of (1), which is a semi-classical field equation, 
contain, by construction, the effects of the quantum fields as represented by 
< Tf^u >■ Therefore, given a solution 

g = 5(0) + hAg (16) 

of (1) (we must distinguish Ag from g, since the former obtains from solv- 
ing the differential equation (1), whereas the latter represents an arbitrary 
fluctuation to be integrated out) we assume that Z contains the action of 
the semi- classical metric as its leading term: 

Z^^-i[9^'^+nA9]^ (17) 

Moreover, since (1) can be obtained from an variational principle |jl^, the 
form of the action / for the semi-classical metric (16) is identical to that 
corresponding to the classical background metric. 
The action is given by |^ 

I = h — I subtract (18) 

where 

Ii = --^ r [ * d''x^R+^ [ d'x^lC. (19) 
167r J2M Jo Svr Js^xS'^ 

The four space metric (16) is given in the Appendix. The boundary at Tq — 
const, is the product x S'^ of the periodically identified Euclidean time 
with the unit two-sphere of area A = Airr'^. The proper length around the 
coordinate is /3* = SvrM. The trace of the boundary extrinsic curvature 
is IC and denotes the induced three-metric. The volume term in (19) is 
sensitive to the trace anomaly < >, as can be seen immediately by taking 
the trace of equation (1). These anomalies have been evaluated exactly for 
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the conformal scalar |1C] and the U(l) gauge boson |11] while an analytic 



approximation has been given for the massless spinor case |12|. They are 

(20) 
(21) 
(22) 



< ^scalar 



< Tj^ ^spinor 



< ^vector 



e w 



As these are all of order 0(e), we can replace ^ — > r^sin^ under the 
integral. Calculation of the volume term is immediate and yields the result 



volume term 



3K 



(23) 



where C G (1,|,— 13) is a spin dependent constant and Wq = When 
integrating (20-22), M is of course held constant, as it depends only on the 
upper endpoint of the integration interval, while w = 2M(ro, T(ro))/r varies 
from 1 to Wo < I- 

The calculation of the boundary contribution to the action is slightly 
more involved. The determinant of the induced three-metric is 



(24) 



The trace of the extrinsic curvature tensor is 
IC 



2 _i /2 1—1/2^1 / N 
-9rr' +T,9rT' ^ log fftf, 

r Z or 



(25) 



Substituting the metric components from Appendix (45,46) into (25) and 
performing the integrations indicated in (19) yields 



boundary term = (87rMr)(l — w) 1 + e{p — riK ^ — w{l — w) ^/i} 



+ 



47rM' 



T.-i 327: ■ 
l + e{p-nK ^ + ^ H ^ < T; >} 



, (26) 
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where use of the equations of motion A(47-48) has been made in order 
to arrive at this final expression. Lastly, the subtraction term is just Ii 
evaluated for a flat four-metric having the same period for the coordinate: 



' subtract 



-(3r, (27) 



where (3 = T ^(r) is the inverse local temperature. Putting all this back 

together, we obtain 

- (87rMr)(l -w)\l + e{p - nK'^ - w{l - w)"V} 



- (47rM2 



l + e{p-nX-i+/i + ^^^<r;>} 



+ Pr, (28) 



where it is understood that this is evaluated at r = ro, and that the functions 
HjP and < >, as well as the constants n,C are spin dependent (see the 
Appendix). 

We are now ready to calculate various thermodynamic state functions of 
interest. In particular, the thermal energy E in the canonical ensemble is 



f diog{z) \ ^/m 

d(3 Ja^KOP 



(29) 

A 



while the free energy F and the entropy S are 

F = -p-^ log(Z) « p-^I, and (30) 

S = /3E + log(Z) ^PE-I. (31) 

In calculating E, there are two possible ways to proceed. On can, for exam- 
ple, compute the derivative in (29) directly. Although straightforward, this 
is somewhat involved due to the fact that the black hole mass M appearing 
in / is a function of the temperature (and cavity radius) and hence, so are e 
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and w = as well. The action must therefore be regarded as a function 
of T (or (3) and r, and thus 

dl\ f dl 2dl 2edl\ fdM\ 

= '— + rir--T7ir: [^] ■ (32) 



dp) A \dM rdw MdeJAKdp) 

As such, a useful ingredient in these calculations is the quantity (^■)a which 
can be derived directly by first differentiating (3 = (using (2)) partially 
with respect to M at fixed radius r, and then inverting the expression so 
obtained, remembering to expand the inverse derivative only to 0(e). We 
find, following this proceedure, that 

/dM\ _ (l-w)V^ j giLli^ 



[dp J A 8^(1 -fw;) 1 +'(l-|u;) 



1 + ^JY^^i^) 



+ (^[f(l--rV + ^<T;>]^ (33) 



where J'{w) = p — nK~^ — ^(1 — w)^^fj,{w). The stress tensor enters here 
because we have employed the semiclassical equations of motion (see the 
Appendix) to eliminate the derivatives of the metric functions. Attention 
should be brought to the pole at w = | (r = 3M), as we might expect this 
to lead to singularities in some thermodynamical quantities at the circular 
photon orbit. A similar pole was found in the expressions for the mass 
perturbations (11). In fact, the specific heat at constant cavity wall area 
Ca = ~P''^{^)a indeed diverges at r = 3M, a result which was already 
known for a vacuum (e = 0) black hole |p. The pole in Ca can be understood 
as arising from the singularity in (33). 

Alternatively, we can make use of the fact that the thermal energy (29) 



is identical to the quasilocal energy, as demonstrated in |14| for gravitational 
systems possessing arbitrary static and spherically symmetric metrics. For 
such metrics, the quasilocal energy is given by jl^ 

E = r - r[g'"' {r)]^ . (34) 
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For the case at hand, we have 



(35) 



with 



= r — r(l — ■u;)2, and 




(36) 



where E^'^^ is the quasilocal 



energy corresponding to a vacuum Schwarzschild 



black hole and /S.E the correction due to the back reaction. 

To get a feeling for the nature of the correction term, we plot the (scaled) 
ratios K-^^ in Fig. 4 (K = 38407r is a constant appearing in the metric 
functions p and n). The necessity of employing a finite spatial boundary 
becomes vividly apparent upon examining the large-r limit of these ratios. 
Inspection of the functions fi{w) in A(53,55,57) shows that fi r'^, whereas 
E^'^^ M. Thus, for any fixed value of e < 1, the contribution to the 
thermal energy coming from the back reaction eventually dominates if r is 
unbounded. This translates into large energy fluctuations and the conse- 
quent instability of the system. The relative corrections for the fermion and 
conformal scalar are plotted in Fig. 4. While the fermion contribution is 
positive definite, the conformal scalar perturbation exhibits a minimum at 
roughly r ~ 2.44M and passes through zero at r ~ 3M. This perturbation 
is negative from 2M < r < 3M and is due to the fact that the effective mass 
function iJ,{w) is negative in the same interval. The renormalized stress ten- 
sors typically violate all the classical energy conditions and the violation of 
the weak energy condition is what is responsable for the negative energy 
correction seen here. Similar behavior is exhibited by the gauge boson, as 
shown in Fig. 5. There, the effect is much larger than in the two former 
cases. (All these corrections have been scaled by the large constant K for 
improved visibility.) 
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In a similar fashion, the entropy (31) of the combined system of black 
hole plus radiation can be split into two contributions: 

S = + eAS, (37) 

where 

= AttM^ (38) 

is the Bekenstein-Hawking entropy of a Schwarzschild black hole and the 
correction term is 



A5 = AttM^ { p{w) - nK-^ + ii{w) + —3— <T^ >] 7^ [l-w^). 



32ttM^ \ USttCM^ 

< T > 

ew^ ' J 3K 

(39) 

The (scaled) fractional correction K{-^^) is plotted in Fig. 6 for the scalar 
and fermion and in Fig. 7 for the gauge boson back reactions. Using the 
explicit forms for the functions appearing in (39) and the indicated compo- 
nent of the stress tensors which are readily available from |10, 11, we 
have 

AS" 2 44 
K{-—-) = -u)-^ + 2u;-2 _^ Qyj-i _ 8 iog(u;) - Ww - Gw"^ + 22w^ , (40) 

for the conformal scalar field, 

,^,A5, 7,4 o „ 1 , , , 180 124 . 92 32, 

^(^) = -[-w-^+4w-^ + 12w-'-mog{w)-—w-—w^ +—w'-—], 

(41) 

for the massless spin- 1/2 fermion, and 

AS* 4 332 
iC( ) = + Aw~^ + 12w'^ - 16 log{w) + 420w - 52w'^ + —w^ - 496, 

(42) 

for the ^7(1) gauge boson. These corrections have the desirable feature that 
AS = at w = 1 (r = 2M). This means one can think of adding layer upon 
layer of entropy, associated with the hole and a given < T/f > beginning at 
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the horizon r = 2M and ending at the cavity wall r = Tq. At w = 1, with no 
"room" for the fields to contribute anything further, one then obtains only 
the Bekenstein-Hawking entropy ^Affh''^, as would be expected. Thus, one 
must regard AS" as arising from both the radiation fields and their effects 
on the gravitational field. These corrections are also of the same order as 
the naive flat space entropy. Again, as is by now familiar, we come across 
an apparent spin dependence in the correction terms. While the fermion 
contribution is positive definite, the fractional entropy corrections coming 
from the conformal scalar and U(l) gauge boson are negative in the ranges 
2M < r ~ 3M and 2M < r ~ 5.5M , respectively. It should be emphasized, 
of course, that the total net system entropy is positive definite for all r > 2M. 
Nevertheless, it is unexpected that the presence of spin-0 and spin-1 fields 
should tend to diminish the entropy in the neighborhood of the horizon. We 
must ascribe this phenomenon to the type of boundary conditions employed. 
When microcanonical (fixed energy) boundary conditions are used, these 
fractional corrections are positive definite in all three spin cases using the 
same set of back reaction solutions |2| . 

4 Black Hole Nucleation 

The free energy in (30) can be used to determine the likelyhood of the 
nucleation of semi-classical black holes from hot flat space. Since the free 
energy is least when the system is in a state of thermodynamical equilibrium, 
the idea is to evaluate F for different phases and then identify the phase 
with the minimum value for F. Since / = (3F and (3 > 0, we can also search 
for minima in the action. Hot flat space is defined as massless quantum 
fields on a flat background geometry. Given the flat space radiation entropy 
Shfs = ^aT'^V and the thermal energy Ehfs = aT^V, where a = 7r^/15/i^, 
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the free energy for hot flat space is 

Fhfs = Ehfs - TShfs = - < (43) 

and is negative. The corresponding action is 

Ihfs = I3FHFS = (44) 

This is to be compared with the action of a semi-classical black hole, Eq.(28). 
This is easiest to do for the high temperature limit where the mass perturba- 
tions for the M2 branch are unimportant (see, e.g.. Fig. 2). In this regime, 
we can approximate M2 as 

to obtain an estimate for /(M2): 

/(M.) = 3.-„^-g+.^(-i^f, (46) 

with b = 96, —32 or 2432 for the conformal scalar, spinor or gauge boson, 
respectively. It is clear that for T ^ 00, the action for hot flat space can 
be made arbitrarily negative (and therefore its free energy) so hot flat space 
should be the dominant phase at high temperature. This conclusion holds as 
well in the limit of zero back reaction [^, but the point here is to determine 
whether the back reaction tends to enlarge or diminish this phase. In fact, 
relative to the vacuum black hole action (or free energy), we see that the 
scalar and U(l) back reactions tend to extend (in temperature) the HFS 
phase, as these correction terms are positive, while the fermion tends to 
diminish this phase. For the low temperature limit, /? — > 00 and Ihfs ~^ 0. 
However, in this regime (a ^ vr in Figures 2 and 3) the mass perturbations 
become important and a higher-order (or non-perturbative) treatment of F 
is probably required to settle the issue. 
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Independently from the temperature dependence, we can explore the r- 
dependence of the action. For the vacuum black hole, the action is negative 
from 2 < r/M ^ 2.25, as can be seen from Fig. 8. This implies F < in the 
same interval. Turning on the back reaction leads to corrections which are 
presented in Fig. 9 for the massless spinor and conformal scalar (scaled up 
by the factor K). The fermion appears to make the action (and hence the 
free energy) more positive relative to the vacuum case, but the scalar induces 
a localized negative component to the action in the range 2 < r/M ^ 3. The 
U(l) gauge boson component is even more striking (Fig. 10), exhibiting a 
substantial negative correction (roughly 60 times deeper than the scalar) 
which extends out to r/M 30. Qualitatively then, we conclude that the 
lowest order back reaction makes the free energy more negative than the 
vacuum case, and thereby increases the possibility for nucleation. 

5 Discussion 

Lowest order solutions of the semi-classical back reaction problem have been 
used in this paper to calculate the 0{ti) contribution to the black hole's 
mass, entropy, thermal energy and free energy in the canonical ensemble. 
We have calculated the seperatc contributions coming from spin 0, ^ and 
1 quantum fields as there are important qualitative distinctions among the 
different spins. Already at lowest order we find evidence for competing 
effects between the half-integer and integer spin field back-reactions in the 
spatial region 2M = r ~ 3M. This is important, because this is also the 
same interval in which the specific heat Ca is positive definite and where, 
consequently, the semi-classical black hole is (locally) stable. The algebraic 
sign differences between the fermion and scalar and gauge boson fractional 
corrections become all the more important if the effect of multiple species 
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are taken into account, such as would arise when couphng the black hole to 
multiplets associated with gauge theories of particle physics. For example, 
the standard model contains 45 fermions, 12 gauge fields and 4 scalars, 
whereas the minimal SU(5) grand unified model contains 45 fermions, 24 
gauge fields and 34 scalar particles. These integers can be used to estimate 
the back reaction induced by multiplet fields, as discussed in ||3|, |^. 

The canonical boundary conditions employed here require a static bound- 
ary held at a fixed temperature. This situation could be dispensed with by 
adopting conditions of fixed pressure and temperature. In this case, the 
appropriate thermodynamic potential is the Gibbs potential G = F — pA. 
This function can be obtained from an appropriate action, as shown in |14], 
by adding the boundary term pA to the canonical action. To make the 
system dynamic, a time-dependent ambient temperature can be introduced 
by embedding the black hole plus radiation into a Friedmann-Robertson- 
Walker cosmology. The external heat reservoir is thus realized in a natural 
way. The boundary seperating the black hole and radiation from the cos- 
mological background becomes dynamic, behaving as a domain wall. The 
embedding could be carried out using the techniques of spacetime surgery. 
At a given time (or equivalently, a given ambient temperature), one excises 
a spherical region from the FRW background and in its place, grafts in one 
of the modified black hole metrics (16) having the same instantaneous ra- 
dius as the excised region and temperature as the FRW spacetime. The 
two spacetimes are joined using the standard junction conditions. If the 
background expansion is not too rapid, one might argue that the system is 
in quasithermal equilibrium with the ambient spacetime. 

Lastly, we comment on the semiclassical back reaction program and 
its possible relation to a (correct) quantum gravity. The solutions of (1) 
employed here contain the 0{h) corrections to the classical background 
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Schwarzschild spacetime. One can also ask for the influence of this met- 
ric on the ambient fields, that is, by computing the radiation stress tensor 
over the new, perturbed, metric (16). This source would now contain terms 
of 0{h^). This source could then be inserted into the right hand side of 
(1) and one could solve for the "back-back-reaction", and so on. In other 
words, (1) is to be solved self-consistently for the equilibrium metric. How- 
ever, while the lowest order solution of (1) is expected to reveal qualitatively 
reliable information, carrying out this iterative program to higher orders us- 
ing only (1) may run the risk of bypassing essential features of a putative 
theory of quantum gravity. This concern is suggested by the fact that the 
semiclassical field equation (1) is derivable from an variational principle, 
which includes the renormalized quantum stress tensor as part of a dynam- 



ical theory involving gravity [15|. In fact, the variation of that (effective) 
action leads to the equation 

G^. + Ag^u + aHj^J + bHj^J = SvrG < T^. >, (47) 

where the tensors Hj^) , H^J are linear combinations of quadratic curvature 
terms arising in the renormalization of the stress tensor. The (renormalized) 
constants A, a, b and G can only be determined by experiment. In partic- 
ular, in order to avoid conflict with observation, it is necessary to assume 
both a and b are very small numerically. To recover Einstein's theory, it is 
necessary to set them identically to zero. The point to be emphasized here 
is that quantum field theory indicates that higher derivative terms are to be 
expected a priori and treatments going beyond lowest order need to address 
this issue. 
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A Appendix 

As shown in Q, the (Euchdean) metric of the perturbed black hole can be 
written as 

(is^ = gtt{w)dt^ + grr{w)dr'^ + r'^{d9'^ + sin^ Odcj?), (48) 

where 

gttiw) = {l-w)[l + e{2{p{w) - riR-^) - wil - w)-^fiiw)}], (49) 

grriw) = {l-wy^il + ewil-wy^niw)], (50) 

and K = 38407r. The two functions /i and p are solutions of the linearized 
semiclassical Einstein equations 

dp 167rM2 



e 



-(l-u;)-i<r;-T*>, (51) 



dw 

dp 327rM2 , ^ ^ 

dw 

These as well as the constant n and the stress tensors actually depend on 
the spin of the quantum field interacting with the hole. Denoting with the 
subscripts S, f, V the conformal scalar, massless spinor and vector boson 
back reactions, respectively, these functions are given by 

Kps = lrdw~^ + 2w~'^ + 6w'^-8log{w)-10w-6w'^ + 22w^-^), (53) 

1 2 40 28 84 

Kps = 2(3^'^ + 4^"^ - 81og(ti;) - - Ww'^ - + y ), (54) 

7,2 o ^ 2 ^ 1 , , 90 62 2 46 , 16, ,,,, 

00 I I o I 

7,2 2 . 1 , / s 200 50 2 52 o 136, , , 
^Pf = ^(^'^ + -81og(t(;)-— w;-yu;'-y?/;^ + — ), (56) 
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Knv = ^w-^ + 2w-'^ + 6w-^ - 8 log(u;) + 210w- 26w'^ + —w^ - 248, (57) 
Kpv = ^w"^ + 4w;-^ - 8 log{w) + + Ww"^ + 4^^ - 32, (58) 
while ns = 12, Uf = —4 and = 304 
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Figure Captions 

Figure 1. The double- valued black hole mass in the canonical ensemble 
plotted as a function of a. 

Figure 2. Mass perturbations Au; for the (heavy) M2 branch arising 
from the spinor, conformal scalar and U(l) gauge field back-reaction. 

Figure 3. Mass perturbations Au; for the (light) M\ branch arising from 
the spinor, conformal scalar and U(l) gauge field back-reaction. 

Figure 4. Black hole fractional thermal energy corrections (scaled up 
by-ftT = 38407r) arising from the fermion and conformal scalar back reaction. 

Figure 5. Black hole fractional thermal energy corrections (scaled up by 
K = 38407r) arising from the gauge field back reaction. 

Figure 6. Fractional entropy corrections (scaled up by K): massless 
spinor and conformal scalar cases. 

Figure 7. Fractional entropy correction due to the gauge boson. 

Figure 8. Action (19) for the vacuum Schwarzschild black hole. 

Figure 9. Back reaction contribution to the black hole action coming 
from the fermion and conformal scalar (scaled by K). 

Figure 10. Back reaction contribution to the black hole action coming 
from the U(l) gauge boson (scaled by K). 
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